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We derive an expression of the core traction contribution to the dislocation elastic energy 
within linear anisotropic elasticity theory using the sextic formalism. With this contribution, 
the elastic energy is a state variable consistent with the work of the Peach-Koehler forces. 
This contribution needs also to be considered when extracting from atomic simulations core 
energies. The core energies thus obtained are real intrinsic dislocation properties: they do not 
depend on the presence and position of other defects. This is illustrated by calculating core 
energies of edge dislocation in bcc iron, where we show that dislocations gliding in {110} 
planes are more stable than the ones gliding in {112} planes. 



Keywords: dislocation theory; anisotropic elasticity; core tractions; clastic energy; core 
energy 

1. Introduction 

One important quantity controlling the physics of dislocations is their elastic en- 
ergy. It is defined as the integral of the elastic energy density over the whole volume 
except a small core region surrounding the dislocation line. This excludes the re- 
gion around the dislocation core where elasticity does not apply because of the too 
high strains. Using Gauss theorem, the elastic energy can be decomposed in two 
contributions: 

• the one corresponding to an integration along the dislocation cut of the work 
necessary to create the dislocation. As it is well known, this contribution varies 
with the logarithm of a characteristic distance of the dislocation microstructure. 

• the contribution arising from the work done by the tractions exerted on the tube 
which isolates the dislocation core. The corresponding contribution to the elastic 
energy is known as the contribution of the core tractions. 

This last contribution is sometimes forgotten. Indeed, the cut contribution is usu- 
ally the dominant one. Moreover, the core traction contribution disappears when 
one tries to define the elastic energy of an isolated infinite straight dislocation 
because of the external cylinder that has to be introduced to prevent the elastic 
energy from diverging. Nevertheless, if one wants the dislocation elastic energy to 
be a state variable, i.e. a variable that only depends on the current state and not 
on the transformation path used to reach this state, both contributions need to be 
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taken into account. Bullough and Foreman pj] already showed that the elastic en- 
ergy of a dislocation loop does not depend on the hypothetical creation mechanism 
only when both contributions are considered. Lothe and Hirth [3] also noticed that 
the definition of the elastic energy of a straight dislocation could not be consistent 
with the work of the Peach-Koehler forces if the contribution of the core tractions 
was not included in the elastic energy. Finally, Gavazza and Barnett 01 showed 
that the dislocation elastic energy part associated with the core tractions leads to 
a contribution to the self- force acting on a dislocation loop. 

Despite its importance, no analytical expressions of the core traction contribu- 
tion to the elastic energy seems to exist in the literature when the elastic anisotropy 
is considered. The anisotropic linear elasticity theory of dislocations has been de- 
veloped in the past sixty years. Eshelby, Read and Shockley [3] were the first to 
express the anisotropic elastic field created by a straight infinite dislocation. Their 
formalism was then enriched by Stroh |^ , leading to what is known as the sextic 
formalism. Latter, it was shown that the elastic field of a dislocation of arbitrary 
shape can be obtained from the fields of straight infinite dislocation^. But no ex- 
pression of the contribution of the core tractions to the elastic energy has been 
obtained until now. The variation of this contribution with a rotation of the dislo- 
cation cut is known [1,0]. It actually corresponds to the angular dependence of the 
interaction energy between two dislocations derived by Stroh [5]. But an expres- 
sion of the absolute value of this contribution is not available. Such an analytical 
expression exists under the assumption that the elastic constants are isotropic 
The few studies that considered the elastic anisotropy as well as the contribution 
of the core tractions l3, 11, l3, 13] calculated this contribution with a numerical 



integration along the surface of the core cylinder. It is the purpose of this article 
to obtain an analytical expression of this contribution in the framework of the 
sextic formalism. Such an expression could be used then in further studies. Poten- 
tial applications are the extraction from atomic simulations of dislocation energy 
properties like core energies [1, d, Q, 15, [H, 17, 111], calculations of dislocation 



loop self energy [l, H, 1^, or computation of the self- force acting on a dislocation 
loop [1] . Such an application is presented in this article where core energies of edge 
dislocation in a-iron are determined. 

In the first section, we define the elastic energy of a straight dislocation so as to 
clearly make appear the contribution of the core tractions. We then use isotropic 
elasticity to illustrate this definition and highlight the importance of this contribu- 
tion for the coherency of the elastic energy definition. The analytical expression of 
the core traction contribution within anisotropic linear elasticity is obtained in the 
following section using the sextic formalism. We finally illustrate the consistency 
of our results by studying edge dislocations in a-iron and calculating their core 
energies with different simulation methods. 



2. Elastic energy of a straight dislocation 

We assume in this article that the elastic field created by a dislocation can be 
reduced to the Volterra solution. Eshelby et al. [1] have indeed shown that a straight 
dislocation in an infinite elastic medium creates in a point defined by its cylindrical 



^For a review, see Ref. 0I- 
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coordinates r and 6 a displacement given by the Laurent series 

oo ^ 

u{r,0) =vln{r)+Uo{d) + J2'^n{0) — . 

n=l 

The two first terms of this series, v In (r)+uo{9) correspond to the Volterra solution: 
this describes the elastic field of the dislocation far enough from its core. The 
remaining terms of the series {n > 1) are the dislocation core field [l^, which 
arises from non-linearities in the crystal elastic behavior and from perturbations 
due to the atomic nature of the core. We do not consider this part of the elastic field 
in our definition of the dislocation elastic energy = Vn > 1) and take only 
the Volterra solution. The elastic field creates nevertheless tractions on the surface 
which isolates the dislocation core. As a consequence, a contribution of the core 
tractions to the dislocation elastic energy exists even when the dislocation core field 
is neglected. As it will be shown below, this contribution needs to be considered 
so as to obtain an unambiguous definition of the dislocation elastic energy. 




Figure 1 . Definition of the contour used to calculate the dislocation elastic energy, (a) Isolated dislocation, 
(b) Isolated dislocation dipole. 

We first consider an isolated dislocation. Two cylinders centered on the disloca- 
tion need to be introduced so as to define the elastic energy (Fig. [1^). The inner 
cylinder of surface Sc isolates the dislocation core. Strains are much too high close 
to the dislocation core to be described by elasticity theory. As a consequence, 
elastic fields are diverging at the origin and one needs to exclude the core region 
from the elastic description. As the elastic energy integrated on an infinite volume 
tends to infinity, one also needs the external cylinder to prevent the elastic energy 
from diverging. The dislocation elastic energy per unit-length of dislocation is thus 
defined in the volume V comprised between both cylinders 

where a and £ are respectively the stress and strain created by the dislocation. 
The use of Gauss theorem allows to transform this volume integral into a surface 
integral: 

E^^o = \ ll^^vUidS,, (1) 
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where 5 is a surface enclosing the volume V and isolating any discontinuity of 
the elastic displacement u generated by the dislocation. A dislocation is by defi- 
nition the frontier of a surface that has been sheared by the Burgers vector b. A 
displacement discontinuity therefore occurs on the dislocation cut which needs to 
be isolated by the surfaces So+ and Sq- (Fig. We would like to stress that 
this cut does not necessary correspond to the dislocation glide plane. The surface 
appearing in the integral [T] is then composed of Sc, Sq- , Soo, and 5*0+. 

The integrals along the surfaces Sc and cancel because the resultant of forces 
located in the core is null, and the normals to the surfaces 5c and 5*00 have an 
opposite orientation [1, @|- As a consequence, the elastic energy of an isolated 
dislocation integrated between the cylinders of radii Tc and Roo is limited to the 
contribution corresponding to the dislocation cut, leading to the the well-known 
result 

j-elas _ 1 1, t-O - - I ^'^ 



Ef(^, = -b.^b,ln^-^j, (2) 

where the tensor K*^ only depends on the elastic constants, and Tc and i?oo are 
respectively the radii of the inner and external cylinders. 

We consider then an isolated dislocation dipole. The first dislocation of Burgers 
vector b is located at the origin and the second one of Burgers vector — b at the 
point defined by its cylindrical coordinates {d,(j)) (Fig. [TJj). The two dislocation 
cuts need to be orientated so that the displacement discontinuity cancel except 
in the surface bounded by both dislocations. This ensures that no displacement 
discontinuity occurs far from the dipole and allows the dipole to be mechanically 
isolated. The elastic energy created by the dipole in the infinite volume is then 



dipole - 2 I L V »i ~^ V j V"' ^ "'i j ^'^1' 

where cr^^^ and cr*^^) are the stresses created by each dislocation, and u^'^'> and tt^^^ 
the corresponding elastic displacements. The integration surface is composed of 
the two cylinders s'^'^ and S'c^'' of radii Tc removing the elastic divergence at the 
dislocation cores, and of the two surfaces Sq- and Sq+ removing the displacement 
discontinuity along the dislocation cut (Fig. Wp)- One does not need to introduce 
an external surface like for the isolated dislocation because the elastic energy inte- 
grated on the infinite volume now converges. 
The integration on both core cylinders leads to the same contribution 

^elas/ ,N _ 1 [[ ^(l),,{l)jc. _ 1 ff ^{2)„(2)^ 



K«i-(0) = - / / a).>uY'dS, = - a .'ul^'dSi. (3) 

2 J 75.(1) ' 2 J Js(/> '1 ' ^ ^ ^ 

The elastic energy of the dislocation dipole is then 

i^dipoie = 2Ef^{ct>) + hKf^b, In . (4) 

We will show in the following that the contribution E^^'^^ of the core tractions 
only depends on the angle (j) defining the azimuthal position of the dislocation 
dipole. It does not depend on the core radius nor on the separation distance 
d. This contribution is not present in the elastic energy of an isolated dislocation 
(Eq. [2]) . This is a consequence of the introduction of an external cylinder to prevent 
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the elastic energy from diverging. But an isolated dislocation cannot exist: another 
defect, like a dislocation with opposite Burgers vector or a surface, is always needed 
to close the dislocation cut and allows mechanical equilibrium. Therefore, one has 
to use the "trick" of the external cylinder so as to define the elastic energy of an 
isolated dislocation. This gives birth to an artefact as the core traction contribution 
E^^^^ then disappears. 



3. Isotropic elastic media 



We first consider that the elastic constants are isotropic. This allows to obtain 
simple expressions of the elastic energy which can be easily manipulated so as to 
illustrate the importance of the core traction contribution. We assume that the 
crystal is oriented in such a way that 63 corresponds to the dislocation axis and ei 
is collinear with the edge component of the Burgers vector (62 = 0). With such an 
orientation, the tensor K*^ appearing in the elastic energy of an isolated dislocation 
(Eq. [2]) or of a dislocation dipole (Eq. U]) writes [IT 



2n 1 




where // is the shear modulus and v Poisson's ratio. 

The core traction contribution has been calculated by Bullough and Foreman [l| 



167t(1 - u) 



1 



1 



2 cos (2(/>) 



One sees that this contribution is null for a pure screw dislocation (61 = 0) , but this 
result does not hold anymore when elastic anisotropy is considered (§11]). For a pure 
edge or a mixed dislocation, if the dislocation cut corresponds to the dislocation 
glide plane, one has then (j) = {) or cj) = n, and one recovers the expression given by 
Hirth and Lothe (Eq. (3-53) in Ref. 0): 



fib 



2 1 



2u 



167T (1 - u) 



2 • 



(5) 



As Poisson's ratio is smaller than 1/2, this shows that the core traction contribution 
reduces the dislocation elastic energy when the dislocation cut corresponds to its 
glide plane. 




Figure 2. Reversible thermodynamic cycle for a dipole of edge dislocations. 
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To understand the importance of the core traction contribution, one can look 
at the variation of energy for a dipole of edge dislocations subjected to the cycle 
sketched in Fig. [2j This cycle has been proposed by Hirth and Lothe (Ref. and 
section 8-3 in Ref. [2Q]). The dipole is first created with both edge dislocations 
lying on the same glide plane and separated by a distance re- The energy required 
to create the dipole is lE^"^"^ + 2Ef^ {(f) = 0), where E^°^^ is the dislocation core 
energy, i.e. the part of energy which cannot be described by linear elasticity. One 
dislocation is then displaced through the cycle 1 — > 2 — > 3. The work performed 
by the Peach-Koehler force during this displacement is {1 — v) /2n. When the 
dislocations are at a distance rc, the dipole is destroyed and one recovers the energy 
2£;core _^ 2Ef'''' ((/) = 7t/2). The variation of energy through the complete cycle is 
then 2^;;:'^'' ((/> = n/2) - 2Ef'^ ((/> = 0) + ^uft^ (i _ j,) = 0. If one had neglected 
the core traction contribution to the dipole elastic energy, the energy variation 
would have been non-null which violates the law of thermodynamics. The proper 
consideration of all energy contributions in the Volterra elastic field thus allows 
a coherent definition of dislocation energetics. Such a coherency does not require 
more complex descriptions which assume a spreading of the dislocation like in the 



Peierls-Nabarro model or in the standard core model proposed by Lothe [2l[ |. 



4. Anisotropic elastic media 
4.1. Sextic formalism 

We consider a dislocation of Burgers vector b located at the origin and we orient 
the crystal so that the axis 63 corresponds to the dislocation line. We assume that 
the angle (p in Fig- H is equal to —n: the dislocation cut corresponds with the half 
plane defined by xi < and X2 = 0. We will generalize in a second stage to a 
different angle (j). 

The elastic displacement and the elastic stress created by the dislocation are 
given to the first order by the Volterra solution. They can be expressed using the 
formalisn{3 developed by Eshelby et al. 0] and extended by Stroh [1, 0]- The dis- 
placement and the stress calculated in a point of Cartesian coordinates (xi, X2, X3) 
are then 



1 ^ 1 

{XI,X2) = -^=F^-4^^alog(a;i +PaX2) 
a=l 



(6) 



a=l 



The sign =p appearing in these equation^ means — for 1 < a < 3 and -|- for 
4 < a < 6. 

The matrices 5^"^ are obtained from the elastic constants Cijki expressed in the 



dislocation reference frame: 



^ijk — Cijkl +PaCijk2- 



^For a more detailled presentation of the sextic formalism, cf. chap. 13 in Ref. [20| of Refs. 7] and [2^ . 
^We use the same sign convention as Hirth and Lothe [20|. Eshelby et al. Q and Stroh [Sj, |3 use the 
opposite sign in Eq. |6] 
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The roots pa are solution of the sextic equations corresponding to the fohowing 
determinant being equal to zero 

\{Br,k+PaBr2k}\ = o, 

and the non-null vectors A'^, associated to each root pa, check the following equa- 
tion 

{B-,+PaBr,,)At = 0. (7) 

In all equations, we use the Einstein summation convention on repeated indexes, 
except for indexes written in Greek letters. When summation on the six different 
roots Pa is required, it will be explicitly written like in Eq. [6l 

The six roots pa are necessary complex. If pa is solution of Eq. [71 its complex 
conjugate pa* is also a solution of this equation. We can therefore assume that the 
six different roots have been sorted in such a way that 

Q{pa) > and pa+3 = Pa*, 1 < a < 3, (8) 

where Q{pa) is the imaginary part of pa- 

The elastic field given by Eq. [6] corresponds to the one of a dislocation with a 
Burgers vector b and a line direction 63 if the constants check the equations: 

a=l 

1 ^ 

2 Bi2kAk^' 
0=1 

We choose the principal determination for the complex logarithm appearing in 
Eq. O The elastic displacement created by the dislocation therefore presents a dis- 
continuity in the half plane of equations xi < and X2 = defining the dislocation 
cut. 

The system of linear equations [9] can be easily solved following Stroh method 
@, [^. To do so, we define a new vector 

= B?2kAk = -Sjifc^fc- (10) 

Pa 

Both definitions are equivalent because of Eq. [71 As the vector A'^ is an eigenvector 
defined by Eq. [71 its norm is not fixed. One can therefore choose it so as to check 
the following normalization condition: 

2AfL^ = l, Va. 

Using the orthogonality properties 0, HI] of the vectors Af and Lf, the 

solution of the system of equations [9] is given by 

Da = -2Lfhi. (11) 



(9) 



0. 



With such definitions, the tensor K appearing in the elastic energy of an isolated 
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dislocation (Eq. [2]) or of a dislocation dipole (Eq. H]) is given by 

a=l 

4.2. The core tractions contribution 

The contribution of the core tractions to the elastic energy (Eq. [3]) is given by 

Ef'^ = -\r [<y^l{rc.0)cos{e)+ai2{r,,e)sm{e)]u^{r,,9)r,de. (13) 
^ Je=-n 

The expression |6] shows that the elastic displacement u{r,9) is the sum of an 
angular function and a term depending only on In(rc). This last term leads to 
a contribution in Eq. [13] which is null because the resultant of forces located in 
the core is null for a pure dislocation 0, 0]. Only the angular part of the elastic 
displacement leads to a contribution in Eq. [T3l One thus obtains 



^elas = _i ^_LAfDa V T—A^D 



c 



27Ti ^ " ^ ' 27ri 

a=l 13=1 



- i^.^.cos (^) + i^Lsin(g) ^ ^.^ ^ 

cos (61) + sm (61) 

This expression already shows that the contribution of the core tractions to the 
elastic energy does not depend on the radius Tc of the core cylinder. Using the 
property [7] of the matrices B^^j^ and the vectors A^, as well as the definition 1101 of 

the vector , one gets 



Ef- = ^ E ±^ E ^^,D^AfL^D,j\p^,p,), (14) 



27ti ^ 27ri 

«=1 13=1 



where the integral (p, a) is defined by 

J\p,q)= r ZP^^^l±^iog(cose + qsme)d0. (15) 
cost'+psint' 

An analytical expression of this integral is obtained in the appendix |X1 Using this 
expression with the properties [S] checked by the roots Pa, one obtains 



Ef-^ = - ^ log (i ± p^) ^ ±—D^ [a^L^ - LfA^) Dp 

a=l (3=1 

^ 3 6 

+ ^ E E (^"^f - ^"^f ) Dp log [Pa - Pp). (16) 
a=l p=A 

We therefore obtained an analytical expression of the core traction contribution 
within the sextic formalism. 
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4.3. Angular dependence of the dipole elastic energy 

We now examine how the dipole elastic energy varies when it is rotated with re- 
spect to the crystallographic axes. The answer could be obtained by rotating the 
dislocations and the elastic constants so as to calculate in the new reference frame 
all quantities needed to express the elastic energy. But a closed-form expression 
of the angular dependence of the elastic energy can be obtained in a fixed refer- 
ence frame. We will show that this expression is consistent with the dislocation 
interaction energy derived by Stroh [11 . 

The angle (j) can now deviate from the value —n considered above. A rotation of 
the dislocation dipole is equivalent to a rotation of the corresponding cut (Fig. [T]). 
Such a rotation does not modify the vectors Af and Lf 0, 113] • As a consequence, 
the tensor K'^ appearing in the dipole elastic energy (Eq. S]) is unchanged. Only 
the contribution E^^^^ of the dislocation core tractions depends on this angle ip. We 
then call AEf^^{cl)) the variation of the dipole elastic energy, with the convention 
A£;eias(_7^) = 0. To calculate this energy variation, it is useful to write the roots 
of the sextic equations in the form pa = ian{ipa). Such a transformation can be 
performed as long as the roots differ from ±i which only happens in degenerate 
cases due to some isotropy of the elastic constant^ 0]. We can now rewrite the 
integral (Eq. [T5]) appearing in the energy contribution of the core tractions (Eq. UT 



tan (ipa — 0) log 



cos (V'/3 - 0) 



COS {tpp) 



The rotation of the cut modify the roots through the relation Pai'P) = tan (tpa — 4') 
(3]. As a consequence, the elastic energy variation is given by 



AEf^{4) = ^ E ±^ E ^^DaAfL^DpAJl^pa^pp), 



2m ^ 2m 

a=\ /3=1 



with 

AJi(Pa,]5/3) 



- tan (V'a — ^) log 
-Pa cos d -|- sin ( 



log 



cos (V>/3 - Q) 
COS (V'/3 - 



log 



COS (V'/^ — 
cos (V'/j) 



„ COS + pa sin 6 

' cos + Pf3 sin ( 



cos cj) + pp sin ( 



log (cos + pp sin 9) 



This integral is calculated in the appendix[Bl With the result of this appendix, one 
obtains 



AEf^{4>) =^ E E ^^DaAfL^^Dp log (cos cp + pp sin , 



o=l P=l 
6 



- E ±^ E -^"^i*^f (^o^^ <t> + Pa sin ( 



2m 

a=l /3=1 



^The degeneracy can be lifted by adding some noise to the elastic constants. 
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Using the orthogonality properties checked by the vectors Af and L", as weU as 
the definition [TT] of Da, one gets 

1 ^ 1 

AEf^^cl)) = -^±— 5iLfL-6,log(cos0 + p„sin</.). (17) 

One recovers an angular dependence which is in agreement with the expression 
of the interaction energy between two dislocations given by Stroh [3], as weh as by 
Cai et aZ. 0.0 

This expression 1171 shows too that the core traction contribution E^^'"^^ is periodic 
of period n, i.e. is insensitive to an inversion of the cut direction. We have indeed 

1 ^ 1 

AEf^'icP + n) = -Y, ±—hLfqh, log (- cos <P - sincP). 

a=l 

Using the property of log {—z), we obtain 

1 ^ 

AEf"^ {<j, + 7t) =AEf^'' ('^) - 4 biL'iL'^bj if < <p < n 

1 ^ 

=AEf^{ct>) + -Y^ biLfL^h, if - Tt < <A < 0. 

o=l 

Finally, the closure properties of the vectors Lf @, 0,111] lead to the desired result: 

AEf'^{(t) + n) = AEf"^ {(/)). 

5. Edge dislocations in a-iron 

So as to illustrate the need to take into account the contribution of core tractions 
to the dislocation elastic energy, we study now edge dislocations in a-iron, and 
determine their core energies E'^°^^. The core energy is the amount of the dislocation 
excess energy that arises from atomic interaction close to the dislocation core that 
cannot be described by linear elasticity. This is therefore a supplementary energy 
contribution that should be added to the elastic energy. One expects that such 
an energy contribution is an intrinsic property of the dislocation: its value should 
only depend on the dislocation and not on the surrounding environment, like other 
dislocations. 

We use the Fe empirical potential developed by Mendelev et al. |23:] in its mod- 
ified version published in Ref. [IJ]. Thanks to the existence of a cut-off for the 
interactions between atoms, two different methods can be used to determine the 
dislocation core energy. One can either work with an isolated dislocation in an infi- 
nite elastic medium (the cluster approach) , or with a dislocation dipole in periodic 
boundary conditions (the dipole approach). In the following section, we will show 



^The expression used by Li et al. based on the work of Cai et al. is different from Eq. 1171 Li 
obtained A_E'^((/>) = ^ ^^_j^Q^biL°' L°'bj"^X{log{cos ip + pasincl>)}, whereas Eg. 1171 can be rewritten 

AE'^{4>) = ^'^—I'i^biLf L°'bj log {cos (f> + Pa sinip)"^ . Both definitions may be equivalent in the case 
of a screw dislocation studied by Li et al, but we could not demonstrate it. 
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that both methods are coherent and lead to the same core energy as long as all 
contributions are considered in the elastic energy. 

a-Fe has a body-centered cubic lattice and the Burgers vector of the most com- 
mon dislocations isb = a/2 (111), where a = 2.8553 A is the lattice vector. Depend- 
ing on the temperature, different slip planes are observed, either {110} or {112} 
planes. It is therefore interesting to study the energetics of both 1/2 (111) {110} and 
1/2(111) {112} edge dislocations, which both can exist. We will therefore determine 
the core energy of these two dislocations. For the 1/2 (111) {110} edge dislocation, 
we use a simulation box with the orientation ei = [lll]/\/3, 62 = [101] /\/2, 
and 63 _= [121] /VE; for the 1/2 (111) {112} edge dislocation ei = [111] /^/3, 
62 = [I2I] /\/6, and 63 = [101] /\/2. In both cases, ei is collinear with the Burgers 
vector and 63 corresponds to the dislocation line. The dislocation glide plane is 
therefore the plane of normal 62 {(j) = in Fig. [1]). 



5.1. The cluster approach 

In the cluster approach, a single dislocation is introduced at the center of a unit 
cell which is periodic only along the dislocation line and with surface in other 
directions. Atoms are displaced according to the anisotropic elastic Volterra dis- 
placement (Eq. E]). Atoms closer from the external surface than the interatomic 
potential cut-off are kept fixed while the positions of inner atoms are relaxed so as 
to minimize the energy calculated with the potential. This therefore simulates an 
isolated dislocation in an infinite elastic medium. A variant of the method consists 



in relaxing atoms at the surface using lattice Green functions [251, |2a, l27|]. This 
may be necessary in ab-initio calculations because of the small size of the unit cell 
that can be simulated. As, in the present work, we use empirical potential, the 
unit cell is large enough so that the Volterra elastic field correctly describes the 
displacements of atoms at the surface and one does not need to use lattice Green 
function to relax them. 





(a) 1/2 (111) {110} 



(b) 1/2 (111) {112} 



Figure 3. Variation with the cylinder radius i?oo of the energy difference AE (Roo) = E (Roo) — 
^biKfjbj In (Roa/rc) showing the convergence to E'=°'^'= for 1/2 (111) {110} and 1/2 (111) {112} edge dislo- 
cations in o-iron (rc = 6 ~ 2.473 A). 



Once inner atoms are relaxed, we calculate the excess energy E (-Roo) given by 
the interatomic potential and contained in a cylinder centered on the dislocation, 
of axis 63 and of radius i?oo- This excess energy is the difference of energy between 
the system with the dislocation and the perfect crystal for the same number of 
atoms. Only the excess energy of atoms contained in the cylinder is considered. 
This can be performed for different values of the radius Roo in the limit of the 



June 29, 2009 19:9 Philosophical Magazine clouet2009'PhilMag 



12 E. Clouet 

Table 1. Parameters defining the energies of both edge dis- 
locations in a-iron. All quantities are given in meV A~ . 



The core energies 


are given 


for a core radius 


= 6. 




^corc 


Ef'"' (0 = 0) 




1/2 (111) {110} 
1/2 (111) {112} 


286 
329.5 


-87 
-89 


370 
378 



maximal radius allowed by the unit cell. According to the section [2] (Eq. [2|), this 
excess energy should vary as 




In this expression, the tensor K*^ is given by Eq. [T2j It is different for the 
1/2 (111) {110} and 1/2 (111) {112} edge dislocations because of elastic anisotropy 
(Tab. [T]). Looking at the variations of the energy difference /S.E {Roa) = E (Roo) — 
^biKfjbj In (Roo/rc), one expects this quantity to be a constant defining E^""^"^. This 
quantity actually varies but rapidly converges to a constant value for an increas- 
ing radius -Roo (Fig. [3]). Variations for small -Roo arises from the dislocation core 
field, which may be important close to the dislocation core but is not taken into 
account in the present approach. This core field needs only to be considered when 
one cannot reach large enough simulation boxes like in ab initio calculations [l^ . 
The convergence of A-E (-Roo) allows us anyway to define a dislocation core energy. 
Values of this core energy are given for both dislocation in Tab. [TJ 



5.2. The dipole approach 





Figure 4. Sketch of the periodic unit cell used to study a dipole of edge dislocations. 

In the dipole approach, full periodic boundary conditions are used. The total 
Burgers vector of the unit cell has to be zero. Therefore a dislocation dipole is 
introduced in the unit cell. The periodic unit cell we used for the present work is 
sketched in Fig. HI Both dislocations composing the dipole share the same glide 
plane and then <j) = 0. 

The energy of the unit cell is minimized by relaxing all atomic positions. Accord- 
ing to section [21 the excess energy of the unit cell should be 

E = 2E^°'^ + 2Ef^ (0 = 0)+ hK^^hj In (^^^ + (18) 
with d = Ix/'i. -Efntlr is elastic interaction energy between the dislocation dipole 
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10000 20000 30000 40000 

Box surface {number of atoms) 




10000 20000 30000 40000 50000 
Box surface (number of atoms) 



(a) 1/2 (111) {110} 



(b) 1/2 (111) {112} 



Figure 5. Variation with the size of the simulation cell of the dislocation core energy deduced from atomic 
simulations of a dislocation dipole in periodic boundary conditions for 1/2 (111) {110} and 1/2 (111) {112} 
edge dislocations in a-iron (rc = b ~ 2.473 A). Symbols correspond to result of the dipole approach for 
different geometries of the unit cell and the horizontal line to the result of the cluster approach. 

contained in the unit cell and all its periodic images. This elastic interaction energy 
is calculated using the method of Cai et al. [8]. One can therefore deduce the core 
energy from the atomic simulations using Eq. [TSl the excess energy E is directly 
given by atomic simulations and all elastic contributions are calculated using a 
code based on anisotropic linear elasticity. 

The obtained core energies are presented in figure [5] for different sizes of the 
unit cells as well as different geometries characterized by the aspect ratio Ix/ly 
All geometries converge with the size of the unit cell to the same limit and the 
converged value of the core energy is in perfect agreement with the value obtained 
in the cluster approach. The table [T] shows that this agreement is possible only 
if one does not forget the contribution of the core tractions to the elastic energy. 
Without this contribution, the core energies would have been about 90 meV 
lower in the dipole approach than in the cluster approach. 



5.3. Stability of an edge dipole 




Figure 6. Variation with the angle <j> of the energy contribution E'^°^'^ + E° {(p) (rc = b). 

Eq. [18] shows that the energy of a dislocation dipole is the sum of a constant 
term (S™'''^ + and a term depending on the distance d between the two 

dislocations composing the dipole. The values obtained from atomic simulations 
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and anisotropic linear elasticity (Tab.[l]) shows that the 1/2 (111) {110} edge dipole 
is more stable than the 1/2 (111) {112} one, whatever the distance d. This result is 
strictly valid only when both dislocations share the same glide plane {(j) = 0). We 
now need to look to the variation of the distant independent term E^°^^ + El^'^^{(j)) 
with the angle (f), i.e. with the orientation of the dipole. 

Fig.[6]shows the variation of E^°^^ + E^^^^ (cf)) with the angle 4> for both dislocation. 
Whatever the orientation of the dipole, this energy contribution is smaller for the 
1/2(111) {110} than for the 1/2 (111) {112} edge dipole. The 1/2 (111) {110} is 
therefore found as the most stable edge dipole. 

6. Conclusions and discussion 

We obtained within the sextic formalism an expression of the core traction con- 
tribution to the elastic energy (Eq. [T6|) . We showed that this expression agrees 
with the angular dependence of the interaction energy between two dislocations 
previously derived by Stroh [H]. This energy contribution is actually important so 
that the dislocation elastic energy can be a state variable consistent with the work 
of the Peach-Koehler forces. 

This contribution to the elastic energy needs also to be considered when one 
wants to extract from atomic simulations core energies which does not depend on 
the simulation conditions. Both the core traction contribution and the core energy 
do not depend on the length scale of the dislocation microstructure. Their physical 
meaning is nevertheless different. The core energy is a dislocation intrinsic property 
which takes into account the fact that atomic interactions cannot be described 
by elasticity close to the dislocation core. On the other hand, the core traction 
contribution is a part of the elastic energy and depends on the positions of other 
dislocations through the angle <j) defining the dislocation cut. 

The application to edge dislocations in iron showed that dislocation core energies 
could be obtained consistently from different simulation approaches, i.e. the cluster 
or the dipole approaches. Both approaches lead to the same core energy when one 
does not forget to take into account the energy contribution of core tractions in the 
dipole approach. This allowed us to conclude that 1/2 (111) {110} edge dislocations 
are more stable than 1/2 (111) {112} ones. 

The contribution of core tractions is also important when defining the elastic 
energy of a dislocation loop as pointed by Bullough and Foreman [1] . For a glissile 
loop, this contribution generally reduces the elastic energy of edge segments: this 
is true at least for an isotropic crystal (Eq. [5]). As a consequence, this would lead to 
loop shapes which are more rounded than with a simple line tension model where 
this contribution is omitted [23]. As this contribution does not depend on the loop 
size, smaller loops should be more rounded than larger ones, in agreement with 
experimental observations [2^ . This energy contribution could therefore explain 
some discrepancies obtained between a simple line tension model and, either ex- 
perimental observations or results deduced from simulations of the loop self 
stress (3o| . 

Gavazza and Barnett [ij] showed that part of the self-force acting on a loop is 
associated with the core traction contribution to the elastic energy. The obtained 
expression of this contribution to the elastic energy (Eq. [TB]) could therefore be 
used in dislocation dynamics simulations to compute the self-force acting on a 
dislocation segment. 

To conclude, we would like to stress that we obtained the expression of the 
core traction contribution to the elastic energy within the framework of the sextic 
formalism. Anisotropic linear elasticity of line defects can also be handled within 
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the integral formalism^, [12] ■ In the integral method, all quantities defining the 
elastic fields are obtained from angular integrals, so that one does not need to 
look for the solutions of the sextic equations. The recasting of our result within 
the integral formalism stills needs to be done. But this may not be possible as the 
expression [16] we obtained does not make appear any known integral. Getting an 
expression of the core traction contribution to the elastic energy within the integral 
formalism may be therefore a challenge. 
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Appendix A. Integral {p, q) 

In this appendix, we calculate the integral appearing in the contribution to the 
dislocation elastic energy of the core tractions: 

7I/ ^ r -P cose + sine , ■ Q^AQ 

J {p,q) = / 7 — log (cosy + gsmfc'jdt/. 

cost^+psmt^ 

An integration by parts leads to 



17T 
-7t 



{Pj q) = ~ [log (cos e + p sin e) log (cos e + q sin e)]l 

, , ^ . ^, — sin6' + gcos6' , „ 

+ / log (cos6' +psm6') -de. 

J_7t cos e + q sin 6* 

We can conclude that 

jHp,q) = -j'{q,p) (Al) 

It is not possible to directly perform the integration, because of the logarithm 
function appearing in the definition of J^{p,q). To circumvent the problem, we 
derive (p, q) with respect to the parameter q 

dJ^{p,q) —p cose + sine sine 

-de. 



dq J_j^ cos 6* + ]? sin 6* cos0 + gsin( 

We thus obtain an integral of a rational function of cos(^) and sin(0). It can be 
integrated using the residues theorem [sit] : 



9J^(p,g) _ 2ni 



dq i + q 



if > and 9(g) > 0, 



1 — q p — q 
27x1 47x1 

H if 9(p) < and 9(g) > 0, 

\+q q-p 

27X1 



1 - q 



if 9(j?) < and 9(g) < 0. 
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Now, we can integrate this result with respect to the parameter q: 

J\p, q)=- 27rilog {i + q) + fi{p) if Q{p) > and Q{q) > 0, 

=27Tilog (i - g) - 47tilog (p - q) + f2{p) if > and "^{q) < 0, 

= - 2m log {i + q) + 4m log {q - p) + fsip) if < and Q{q) > 0, 

=27Tilog (i - g) + U{p) if < and 9(g) < 0, 

where /i, /2, /s and /4 are four functions depending solely on the parameter p. 
They can be determined by using the property lAll and the result J^(i, — i) = 27T^. 
This leads to the final result: 

{p, q) =2m log (i + p) - 27ti log (i + q) if > and Q{q) > 0, 
=27ri log (i + p) + 27ti log (i — q) 

- 47tilog {p - q) if > and "^{q) < 0, 
= — 27ti log (i — p) — 27Ti log (i + q) 

+ 47Ti log (g - p) if < and 9(0') > 0, 

= - 27Ti log (i - p) + 27Ti log (i - q) if 9(p) < and Q{q) < 0. 



Appendix B. Integral AJ^{p,q) 

The angular dependence of the dipole elastic energy makes appear the following 
integral 



AJi(p,g) 



" —p cos 6 + sin ( 



cos + p sin 

d^. 



, costy + gsin„ . 
log \ — 7 — log (cos + g sm >i) 



cos + (7 sm ( 

We use the same integration method as for {p, q). A derivation with respect to 
the parameter q leads to 

9 ^ , . —pcosO + sinO —s\ii(f) 
—AJAp,q)= / ■ -aO. 

oq ^ J _^ cosa+psm.6 cos<p + qsmcp 

Using the residues theorem, we obtain 

^A4ip,q) =2m if 9(p) > 0, 

oq ^ cos (p + q sm (p 

= _ 2m ^ if 9(p) < 0. 

cos (p + q sin cp 

We then integrate with respect to the parameter q and we use the property 

Aj]{p,q) = —AJl{q,p) deduced from an integration by parts of the initial in- 
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tegral. This leads to the final result 



(P) q) =27ti log [cos + g sin (f)] 



— 2ni log [cos ^ + p sin cf)] 
=2m log [cos + g sin 4)] 



if Q{p)>0 and Q{q) 



>0, 



+ 27ti log [cos ^ + sin 4>] 
= — 2ni log [cos (p + q sin (/>] 



if > and Q{q) 



<0, 



— 27Ti log [cos (f> + p sin (/>] 
= — 2m log [cos 4> + q sin 



if < and Q{q) 



>0, 



+ 27ti log [cos (p + p sin ^] 



if Q{p) < and ^(g) 



< 0. 
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